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ABSTRACT 

This paper develops theoretical models to predict 
the radiation of low-frequency and impul s i ve sound 
from horizontal-axis v/ind turbines due to three 
sources: (l) steady blade loads; (2) unsteady 

blade toads due to operation in a ground shear; 

(3) unsteady loads felt by the blades as they 
cross the tower v^ake. These models are then used 
to "predict the acoustic output of MOD- 1, the large 
wind turbine operated near Boone, N.C. Predicted 
acoustic time signals are compared to those actual- 
ly measured near MOD-1; good agreement is obtained. 

^ INTRODUCTION 

The possibility of satisfying some fraction of our 
energy heeds with wind power has recently received 
serious cons I derat ion. As a result, a number of 
exper imental vertical- and horizontal -axi s wind 
turbines have been constructed under DOE sponsor- 
ship and are presently being tested. In this 
paper we present aerodynamic and acoustic models 
developed or adapted to predict sound radiation 
from hor I zontal -axi s wind turbines due to the 
following three acoustic sources: (l) steady 

blade loads, (2) unsteady blade loads due to wind 
shear, and (3) unsteady loads acting on the blades 
as the latter pass through the tower wake. 

Below, we describe the mathematical models used 
to predict noise due to these sources and apply 
these models to calculate the acoustic output of 
the DOE/NASA MOD- I wind turbine (for which some 
tower wake data and sound measurements are avail- 
able) under typical operating conditions: a free 

stream c f 35 MPH and rotor speeds of 23 and 35 RPM. 
Although possibly an important cause of impulsive 
low-frequency sound radiation under certain oper- 
ating conditions, cross flow into the rotor was 
not studied here due to difficulties in aerodyna- 
mic modeling similar to those met in non-rigid 
wake analyses for helicopter rotors (ref. l). 

Study of this mechanism should be part of future 
theoretical and experimental research on wind 
turbine noise. In the present preliminary study 
we are niaTniy concerned with modeling the rotor 
aerodynamics in order to determine the strength 
and frequency content of sources of sound on the 
rotor disk; we do not consider terrain and atmos- 
pheric effects on the propagation of such sound 
once it has been created. So, the models devel- 
oped here for mechanisms ( 1 ) - (3) , above , do not 
take Into account the presence of the ground plane 
Qj- the effects which the wind profile or an atmos- 
pheric temperature gradient will have on acoustic 
propagation. 

The acoustic models used to predict sound due to 
(1), steady blade forces and (2), unsteady blade 
forces to ground shear are based on, respectively: 
the classical Gutln propeller noise theory (ref . 2); 
and Its generalization for the case of unsteady 
blade sources (ref. 3). In both studies we used 


lifting-line theory coupled with a vortex-lattice- 
wake model to calculate loads at chosen radial 
stations for each blade. The wake was assumed to 
be semi-rigid in that it was allowed to move with 
the local wind velocity (ref. 4). We did not 
assume compactness of acoustic sources on the 
rotor disk which would have simplified the analy- 
sis somewhat. 

To Investigate the acoustic effect of (3), blades 
cutting through the tower wake, we applied the 
blade-slap theory of ref. 5- Here, the unsteady 
airloads acting on a blade passing through the 
mean profile of the tower wake were obtained 
using Filotas’ linear, unsteady aerodynamic theory 
(ref. 6). The wake data used in the calculations 
Is that given in ref. 7 for a l/40th-scale model 
of MOD-1, Predictions of sound are in good agree- 
ment, both quantitatively and qualitatively, with 
field measurements. 

FORMULATION 

A. Sound due to Steady Blade Loads (Gut in noise) 

Even if the tower and rotor wakes did not interact 
with the turbine blades, and ground shear and 
cross flows were absent, sound would radiate from 
a wind turbine due to the rotary motion of steady 
blade loads. The model used here to study the 
contribution of this source mechanism to the total 
far-field signal is based on an adaptation of the 
classical Gutin propeller noise theory (ref. 2); 
propagation of sound was assumed to take place in 
a stationary medium and the acoustic effect of the 
presence of a ground plane was neglected. The 
source representation for each blade did not take 
advantage of acoustic compactness In either the 
chbrdwise or radial directions. Instead, we con- 
structed a distributed source region for the rotor 
plane (zero coning angle) by using values of blade 
thrust and in-plane forces for six radial stations. 
Furthermore , at each station these values of force 
were distributed along the local chord assuming a 
flat-plate loading distribution. 

The actual calculations of blade forces were per- 
formed using the lifting-line aerodynamic numeri- 
cal TTiodel described In detail In ref. The model 
uses curved grids or lattices of trailing and shed 
vortices to represent the wake behind a turbine 
rotor In a nonuniform incident free stream. Each 
lattice has its origin at the trailing edge of a 
rotor blade and is allowed to change shape as each 
of its grid points Is convected with the local 
value of the nonuniform free stream velocity; such 
a wake is “semi -r I gi in current aerodynamic ter- 
minology. The blade loads are calculated using an 
iterative procedure. The code, named US-21 by Its 
original developers (ref. 4), outputs values of 
blade forces around the rotor for an input ground 
shear; it can calculate steady blade loads for the 
special case of a uniform rotor inflow. Although 
based on potential theory, US-21 incorporates some 


401 



viscous effects in the aeroacoustic model by com- 
puting values of local skin-friction drag from 
knowledge of the calculated potential flow field 
and the drag coefficients for each local blade 
section. Figure 1 shows the coordinate system 
used in the acoustic model for the steady loads; 
since the Gutin sound field is axi symmetr i c , only 
the angle \p is needed to describe the predicted 
directivity patterns. 


In the classical theory of Gutin noise (ref. 2) 
the pressure level P at the s^^ harnx>nic of the 
blade-passage frequency is given by the following 


expression: 


PjCr.i/i) =-^ e ° (-i) 




s^cosjL f 


dr' I(r')Ajr')J 


sB 


/sBQr* sin \Ij\ I 

* ^ c / I 


[1) 


where r,i/^ define the observer's position in the 
axisymmetric (in y) Gutin sound field; B stands 
for the number of blades, for the speed of 
sound, and and R for the rotor's speed and 
radius, respectively. T(r') and I(r') denote the 
values in pounds per linear foot of steady thrust 
and in-plane forces, respectively, computed by 
US-21 at six radial stations for each blade; r' is 
the distance separating each station from the hub. 
We chose stations located at 20, 30, 40, 60 , 75, 

85 , and 95^ of the rotor radius. At each, we 
distributed the values of the calculated forces 
along the local ch ord c(r ') us in 
loading shape Ap~/(x+x )/(x-x. 
coefficients of the load shape f 
designated A (r‘) in (1). With y (r‘) indicating 
the position of the trailing edge of each blade 
section for each radial station relative to that 
of r'=0.2R, we have, for the thrusts (the in-plane 
forces were treated similarly) 


g a flat-plate 
) ; the Fourier 
unctions are 


Ap{r' ,y) = 


T(r') ^ 

TTc ( r^ /2 V y-y^ ( r ' ) -c ( r * ) / r ' 

for T^(r')<Y<’('^(r') + 


0 otherwise 


and so, by Fourier’s theorem, 

ao 

Ap(r' ,Yi t)“T( r ' ) J (r ' ) exp[2Tri s (Y-f2t) /F] 
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(3) 
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where t is time. 


A Fast-Fourier Transform routine calculated the 
A^(r‘) coefficients in (3) and the time signal 
from (4). A sample size of 2048 points was used 
to represent the load-shape functions for every 
radial station. The same number was used to obtain 
the time signal from the calculated harmonics. 

B. Low-Frequency Sound due to Operation in a 
Ground Shear 

Since the turbine operates in the earth’s boundary 
layer, the actual free stream can be far from uni- 
form and the associated acoustic field far from 
axisymmetric. Low-frequency sound then has two 
origins: (l) the steady loads in rotary motion, 

as for uniform flow; and (2) the unsteadiness due 
to the variation of blade forces around the rotor. 


As in the Gutin noise study, the acoustic model we 
used to investigate the effect of wind shear does 
not take advantage of compactness of acoustic 
sources on the rotor. The distributed source re- 
presentation was constructed as previously de- 
scribed, that is, blade loads (thrust and in-plane) 
were calculated at six radial stations and at each 
such station their local values were distributed 
across the local chord using the flat-plate loading 
shape. The acoustic theory used is the standard 
generalization of the Gutin model to allow blade 
forces to vary around the rotor (ref . 3) . We 
assumed that sound propagates in an infinite sta- 
tionary medium with no ground plane. Program 
US-21 calculated the blade loads for the turbine 
operating cases of interest. 

We considered an example of a free stream with a 
shear linear in profile constructed from the sum 
of a uniform 35 MPH free stream and a linear func- 
tion with a magnitude of 8 MPH at the highest and 
lowest points on the rotor. Such a model for 
shear, besides approximating well most free-stream 
profiles, enabled us to compare the predicted 
acoustic fields to those obtained here using 
Gutin’s theory because, for the small blade-pitch 
cases considered here, the pressure far field for 
linear ground shear contains approximately the 
same amount of Gutin noise as its uniform-flow 
counterpart. This allowed us to isolate the 
acoustic effect of unsteadiness due to the non- 
uniformity of the incident rotor inflow. 


The standard generalization of Gutin's theory for 
nonuniform incident inflows contains the following 
expression for the harnfK>nics of the acoustic signal 
in the far field (ref. 3): 


tsBf^r . ^ 

— - — +i5By CO 
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Here T is 2 tt/B; for the two-bladed MOD- 1 , r = 7T. 

Having calculated the harmonics by means of (1), 
we obtain the acoustic signal p using 
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where, since the flow into the rotor is not uni- 
form, the thrust and in-plane forces T and I be- 
come functions also of l^he in-plane rotor angle. 
We then write, for the thrust, 


Ap(r' , y)= 


T(r',Y) , / 

TTc ( r ^ /2 V ) “c(r ' )/r ‘ 


for Y^(r ' )<Y<Y^(r' ) 



0 

or otherwise 

so that in place of (3) we now have 
- - “ ■ : ^ 27TlYq ^ 

Ap(r'jX,t)=f B^(r')e I A^(r') 
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where the substitution T(r' ,y)= I B (r’)6 

q=-oo H ^ 

has been made. Here, T^-2'n, and the 8 (r*)'s are 
the Fourier coefficients of T( r ' ,v) > ca 1 cul ated 
numerically for each radial station using a Fast- 
Four ter Transform routine. Equations similar to 
(6) and (7) were also used to express the in-plane 
forces I(r*,Y)* T(r*,Y) and I(r’,Y) changed 

around the rotor only very gradually for the case 
of shear st udie d here; so we found that 

they and their transforms could be adequately re- 
presented by 2k sample points. 

The time signal Is still given by eq. (4). 

C. Impulsive Sound due to the Passage of the 
Blades through the Tower Wake 

In order to investigate sound radiation due to the 
tower wake we applied a modification of the heli- 
copter-blade slap theory of ref. 5- Briefly, the 
present model performs the following sequence of 
calculations: (1) the measured mean tower-wake velo- 
city profile is expressed as a sum of upwash gusts 
by the Fourier transform theorem; (2) the unsteady 
force acting on a blade passing through a single 
gust is determined using Fi lotas' aerodynamic 
theory; (3) the unsteady force acting on a blade 
as it encounters the tower wake is obtained by 
superpos i t ion of solutions from (2); (4) spanwise 
superposition of solutions of (3) is applied to 
rTK)del the three-dimensional effects of blade rota- 
tion and finiteness in span; (5) the wave equation 
in an Infinite stationary medium (with no ground 
plane) Is solved with the aerodynamic forces as 
boundary condition. 

In step (1) above, we used the wake profile appro- 
xifTsately at the 7k% rotor radius station to repre** 
sent the profiles at all radial positions. The 
actual tov/er wake data for the numerical calcula- 


tions was taken from ref. 7, where the 74^ radial 
position corresponds to a height of 23.2" from 
the base of a l/40th-scale model of MOD-i . The mean 
wake velocity measurements were made at points 
downstream of the rotor which were at a distance 
from the tower-center line of 13.13“ (43*77 feet 
in scale); we assumed these did not differ greatly 
from those on the actual rotor plane, which for 
the model Is at a distance of 9.10" (or 30.33 feet 
In scale) from the tower-center line. Two sets of 
calculations were performed, one using a rota- 
tional speed of 35 RPM and the other 23 RPM, In 
each we Investigated wakes for free streams orien- 
ted at 0°, 20°, and 45° with respect to the tower 
structure; for convenience, computer generated 
reconstructions of these velocity profiles norma- 
lized by the free stream are presented in figs. 
2a-c. In each, the abscissa, which represents the 
distance across the tower wake normal to both the 
tower and the free stream, has been scaled to 
actual MOD-1 dimensions. 

To calculate the local unsteady load acting on a 
blade section passing through the tower wake, we 
use Fi lotas* linear two-dimensional aerodynamic 
theory (ref. 6), The model is adapted from ref. 5, 
where it was used to predict noise due to blade- 
vortex interaction in helicopter rotors. Its 
application here i^ justified by the following tv^ 
conditions on turbine operation and tower wake 
characteristic, respectively, which generally are 
met by wind turbines: (l) small advance ratio 

(small ratio of rotor inflow velocity to rotor tip 
speed); and (2) small sector of rotor azimuth 
occupied by the tower wake. 

We assumed the wake was two-dimensional, that is, 
that the same wake profile existed at all posi- 
tions along the tower height. Als o, we used a 
local two-dimensional model for the blade unsteady 
aerodynamics. Later, we incorporated some three- 
dimensional effects of blade rotation and flnlte- 
ness of blade length Into both aerodynamic and 
acoustic models. 

Let w(y) denote the upwash felt by a Blade passing 
through the tower wake; the variable y measures 
the distance normal to the radial direction on the 
rotor plane. By Fourier's theorem, we express 
w(y) as a sum of sinusoidal gusts of varying 
ampl i tude W(s) 

oo 

w(y) = [ ds W(s) e' 

-> (8) 

— oo 

The unsteady load L(t,s) acting on the blade due 
to the interaction with a single gust is obtained 
applying Fi lotas* aerodynamic theory: 


L(t,s) = 


TFp (QR) 

/ 1+2TT5 


c exp 


-7T^/(2+47rs) 


W(s)0 


-1st 


( 9 ) 


where p is the background air density, 0 and R the 
rotor speed and rotor radius, respectively, and c 
the blade chord. We shall refer to the term in 
brackets in (9) as L (s) . From (9), and using 
superposition, we ob?ain the unsteady blade loading 
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L(t) due to the passage through the tower wake 


L(t) 


00 

» I ds L (t ,s) 

— oo 


( 10 ) 


Here, r Is the distance between the observer in 
the far field and the turbine hub, and z is the 
projection of this distance on a plane normal to 
that of the rotor, M is QR/c , the tip Mach 
number. ° 


The expression above (valid for a locally two- 
dimensional flow) has been derived for a blade of 
infinite length traveling rect i li nearly. The 
three-dimensional aerodynamic and acoustic effects 
of blade rotation and finiteness of blade length 
are modeled as follows: the acoustic pressure p 

satisfies the three-dimensional stationary medium 
wave equation; the unsteady loads on the blade 
passing through the tower wake are taken to be 
equal to the two-dimensional sectional load given 
in (10) with the local blade velocity (and thus 
the magnitude of the loading) varying linearly 
from hub to tip. Given the assumption we have 
made here about the tower wake being the same for 
all points along the vertical (from blade tip to 
hub), such a model for the actual three-dimensional 
loading on a rotating blade, which sees a relative 
free stream linearly increasing in the tip direc- 
tion, should be reasonably accurate except at a 
small region near the tip; here it overestimates 
the actual dipo le strength since the latter must 
vanish as /R-r ‘ . We solve the three-dimensional 
wave equation subject to the boundary condition 


Ap(r‘ ,y) = 


t L(t)6{f (y-t)} 


for 0 < r ' S R 


0 


for r' > R 


( 11 ) 


where, as before, r’ stands for the radial distance 
from the hub measured on the rotor plane. All 
spatial variables have been normalized by the blade 
semichard c/2, and time by (c/2)/^R. 


The solution of the above boundary value problem 
has the following form in the far field: 


.(r.f,^,-!p^;t)= I ds s L^(s)D(s;f-,^,^) 


• expj^is (Mr-t) /( 1 -M^) 


( 12 ) 


wi th 




I . 


exp[- isRM(r'/r)/(1 - tt^) ] -1 
R [M(r'/r)s/(l -H^)]^ 


i exp [-isRM(r'/r)/(1 - M^) ] 
s[M(r-/r)/(l -Mp] 

(13) 


The directivity factor D(s;z/r,y/r,r ' /r) in (13) 
reaches a maximum In magnitude when rVr = 0, that 
is, at points in the acoustic field which coincide 
with the rotor axis. Naturally, for a given r 
there are two such points — one upstream of the 
turbine and one downstream. The fourier trans- 
form W(s) , of the tower-wake upwash on the blade, 
was computed numerically with 20A8 sample points 
using a Fast-Fourier Transform routine; so was 
the acoustic pulse in the far field as given by 
( 12 ). 


DISCUSSION OF RESULTS 

A, Predictions of Gutin Noise 

Figure 3 shows predicted directivity patterns at 
one kilometer from MOD-1 for the first three 
harmonics of the acoustic pressure for an opera- 
ting condition of 35"RPM rotor speed and a 35"MPH 
free stream; the sound pressure level had a 
maximum value of 56 dB In the direction i|;=65‘’v 
Figure k corresponds to the 35-MPH free stream, 
23 -RPM condition, for which the maximum level was 
50 dB. A value of zero pitch at the blade tips 
was used in both sets of calculations. 

For both uniform-flow operating conditions studied 
here we found that the predicted spectrum had the 
expected feature of containing only a few harmo- 
nics of significant level; therefore, the time 
signals (not shown here) were essentially pure 
sine waves at the blade passage frequency. Also 
as expected, the predicted Gutin directivity 
patterns indicated zero on-axis sound and a pre- 
ference for radiation in the upstream direction 
of the turbine. Finally, we also confirmed that 
the case of lower rotational speed (23 RPM) 
radiates less overall Gutin noise than that of 
35 RPM. 

B. Predictions of Sound due to Ground Shear 

Because the ground shear produces unsteady blade 
loads that vary with position around the rotor 
disk, the sound field is a function of both the 
angle and angle in the rotor plane y* Normally, 
for an observer on the ground far from the rotor, 

Y will be essentially 90°. However, if the obser- 
ver Is at a position in the field which Is close 
to the rotor or much higher or lower than the 
horizontal plane of the turbine's base, acoustic 
predictions at other values of y may be of Interest. 
Our numerical results indicated however that the 
dependence on y very weak. 

Figure 5 shows the predicted directivity pattern 
at one kilometer from MID-I for the first three 
harmonics of the acoustic pressure for y=90°; 
the maximum sound pressure level was found to be 
about 58 dB in the direction of (|^-65°. The oper- 
ating condition for the turbine was a 35“RPM rotor 
speed, zero blade tip pitch, and the previously 
described linear-shear free stream. Figure 6 
shows similar predicted results for the operating 
condition of 23 RPM, zero blade tip pitch, and the 
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same ground-shear profile. The maximum sound 
pressure level was found to be about 52 dB for 
this case- 
in conclusion, in contrast to the acoustic direc- 
tivity patterns for Gut in noise, those for ground 
shear indicated non-zero values of on-axi s sound 
levels, due mostly to the higher harmonics. This 
relative growth of the higher part of the spectra 
is apparent also at the field points where the 
ca IcuTa ted "sound pressure levels had maximum 
values In that the shapes of the acoustic signals 
there (not pj^e^ented here) had slightly steeper 
slopes than those for Gut in noise. Accordingly, 
the maximum sound levels were higher than those 
obtained for uniform rotor inflow; the difference 
came from the added contribution to the sound 
field of the unsteadiness of blade forces, an 
acoustic source which of course is absent when the 
Incident rotor inflow is uniform. 

C. Predictions of Sound due to Tower-Wake/Blade 
1 nteract ion 

Figures 73-c for 35 show predicted acoustic 
signals at one kilometer from I^D-I for the wind- 
tower angles of 0°, 20°, and for which we ob- 

tained predicted sound levels of 85, 83, and 8^ dB 
respectively, at the field point of maximum acous- 
t Tc ~pre s sure (0° azimuth) . These values are in 
good agreement with the amplitude of the measured 
signal shown in Fig. 8 (ref. 8). The signal In 
Fig. 8 was recorded approximately one kilometer 
from MOD-1 and at a position in the field 2A° in 
azimuth off the turbine's axis (see Fig. 10 below); 
by definition the direction of the wind corresponds 
to a 0° azimuth. The predicted time signals for 
the 23 ”HPH case were qualitatively similar to those 
shown here for 35 8PM but with the substantially 
lower values for maximum amplitude of 7^, 72, and 
73 ciB for wind-tower angles 0°, 20°, and A5°, 
resp^ec t i vef y . Figure 9 shows the predicted sets 
of spectra of these computed acoustic signals for 
the two cases of 35 and 23 RPM. Figure 10 shows 
the predicted simple acoustic dipole directivity 
pattern at one ki Ipmeter for the case of 35 8PM 
and wind- tower angle of 0°; the position in the 
field corresponding to that where the acoustic 
measurement in Fig. 1 was taken is indicated, 

CONCLUSIONS 

Based on the predictions of the models, we conclude 
that neither steady blade loads nor loads due to 
operation in ground shear contribute substantially 
to the acoustic signal from a wind turbine such as 
MOD-1. Hov^ever, comparison of the theoretically 
predicted signal for noise fr<xn interaction with the 
mean wake ( F i g, 7s) fora 35“8PM rotational speed, 
35-MPH wind speed, 0° wind-tower angle, and the 
measured one shown in Fig. 8 indicates a close 
resemblance, both quantitatively and qualitatively. 
The condor ton to be drawn from this that high 
level impulsive sound could radiate from MOD-1 due 
to the interaction of its blades and tower wake. 

The results for the 23‘'8PM case suggest that a 
possible way to reduce a turbine's acoustic output 
is to operate at a lower rotor speed. Instead, 
dramatic noise reduction could be achieved without 
compromising performance fay changing the structural 
configuration of towers to have smoother wakes. 


To accomplish this, the relationship of complicated 
tower geometries and their wakes will have to be 
Investigated experimentally in follow-on studies. 
Earlier we had pointed out that cross flows into 
the rotor might have significant acoustic effects. 
Investigation of this possibility should also form 
part of future experimental and theoretical 
research. 
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NOMENCLATURE 


A Fourier components of the flat-plate load 
^ shape functions, Eq. (l) 

B number of blades - two for MOD-1 

Fourier coefficients of I, T, respectively, 
^ given in Eq. (5) 

C local blade chord 

c sound speed 

o 

D directivity factor, Eq. (13) 

I magnitude of computed in-plane blade force 

J Bessel function of order n 

_n 

L unsteady lift acting on a blade passing 

through a gust, Eq. (9) 
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magnitude of unsteady lift due to blade/gust 
interaction given by Filotas* theory, term in 
brackets in Eq. (9) 

blade tip Mach number 
acoustic pressure 

harmonics of the acoustic pressure signal p 

distance of observer in the far field from 
turbine hub 

radial distance measured on the plane of the 
rotor disk 

rotor radius 

harmonic counter 

magnitude of computed thrust blade force 
time 

Fourier transform of w 

downwash felt by the blades as they cross the 
tower wake 

distance measured across the tower wake normal 
to the free stream direction and tower center- 
1 i ne 

distance normal to rotor disk 


Coordinate 
and ground- 

system used in the Gut in 
-shear model s . 

c 



Figure 3: D 

t 
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attern a 
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Figure 9: Spectra of predicted acoustic signals 

at 1 Km for 35 RPM; wind-tower 

angles 0®, 20°, 45**. 



Figure 10; Predicted acoustic directivity pattern 
for tower -wake model for wind -tower 
angle of 0°. Location of measurement 
is Indicated in the acoustic field. 
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QUESTIONS AND ANSWERS 



W.L. Harris 



From: F.W. Perkins 

Q; Will any internal structural changes to the blades reduce impulsive noise signifi- 
cant ly? 

A. /Vi?, S't'PliQ ohCLtlQ & S 'l-Tl ~t}l& h T CXd& S j WVtl^'l-Yl ^'LTf\X "td't'UOyiS of ^dZ-Tl'tcziflt'l^Q T* 6 Z ‘Vdh T 6 

blades y will not influence the tower wake, A small reduction in the amplitude of 
impulsive noise may be observed if a ^'softer^^ blade is used. 
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